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Abstract. The orbital magnetic moment of a Co adatom on a Pt(lll) surface is 
calculated in good agreement with experimental data making use of the LSDA+U 
method. It is shown that both electron correlation induced orbital polarization 
and structural relaxation play essential roles in orbital moment formation. The 
microscopic origins of the orbital moment enhancement are discussed. 
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1. Introduction. 

According to Hunds rules, gas-phase transition metal atoms possess large spin Ms and 
orbital Ml moments mediated by intra-atomic Coulomb interactions. In a solid, where 
electron derealization and crystal field effects compete with Coulomb interactions, 
there is a substantial decrease in M$ and partial or total quenching of Ml ■ 

Recent X-ray magnetic circular dichroism (XMCD) measurements [I] report Ms 
and Ml of Co-adatom and small Co clusters on Pt(lll) in an ultra-high-vacuum. 
The M L = 1.1 ± 0.1 fi B and M s + 7M D = 1.8 ± 0.1 [i B (where M D is a spin 
dipole moment) were evaluated from XMCD spectra using the conventional sum 
rules [3J. The number of holes in the Co atom 3d- manifold nd=2.4 was taken from 
local spin-density calculations (LSD A). The XMCD experiments are complimented by 
spin-polarized-relativistic Korringa- Kohn-Rostocker (KKR) Greens function LSDA 
theoretical calculations, and Ms — 2.14 /is and A/l=0.60 hb for Co site were obtained. 
No structural relaxation of the Co atom position over the Pt surface was considered, 
and the atomic-sphere approximation was employed. 

The authors of Ref.p] assumed that the Ml discrepancy between the KKR-LSDA 
theory and the XMCD experiments originates from the lack of orbital polarization 
(OP) in LSDA. They used the well known orbital polarization correction of Brooks [3], 
adding to the LSDA total energy functional an ad-hoc term \BrM\ with the LSDA 
calculated Racah parameter Br. This form of OP correction was widely used in the 
past to improve upon Ml in the bulk transition d- and /-metals where sometimes it 
works well [4]. 

For the Co-atom on Pt surface, Gambardella et al. found out that the Brooks 
OP yields the Ml which substantially exceeds the experimental value. They had to 
reduce the LSDA calculated Br by 50 % in order to obtain a Ml=1.50 /j,r comparable 
with the experimental XMCD data. Ii is assumed in pQ that the reduction of Br 
compensates for a lack of structural relaxation. 

In this work we explore another avenue for the orbital polarization correction to 
LSDA which is based on the correlated band theory LSDA+U method 5j. It consists of 
LSDA augmented by a correcting energy of a multiband Hubbard type and a "double- 
counting" subtraction term which accounts approximately for an electron-electron 
interaction energy already included in the LSDA. Minimization of the LSDA+U 
functional generates not only the ground state total energy, but also one-electron band 
structure energies and spin-orbital states. The basic difference between LSDA+U 
method and the LSDA is its explicit dependence on on-site spin and orbitally resolved 
occupation matrices. The LSDA+U method creates in addition to the spin-only 
dependent LSDA potential, the spin and orbitally dependent on-site "+[/" potential 
which gives OP beyond that given by the LSDA (where it comes from the spin-orbit 
coupling only). 

It was shown by Solovyev et al. [6] that LSDA+U produces the correct OP 
for insulating 3<i-oxides. Recent parameter-free GW calculations for transition metal 
based materials [7] produce OP which is very similar to LSDA+U results with the 
appropriate choice of Coulomb-C/ [8] . 

2. Results and Discussion. 

We performed supercell calculations to model a Co adatom at a Pt(lll) surface. The 
supercell consists of three Pt(lll) layers with doubled (p(2x2)) 2-dimensional unit 
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cell, and the Co atom on the top taken in the fee position (see Fig. 1). The vacuum is 
modeled by the equivalent of two empty Pt layers. All in-plane inter-atomic distances 
are adopted to be those of pure Pt. The distance between the Co atom and the Pt 
surface was varied in the calculations. We note that while the chosen supercell is 
quite small, it provides separation of Co atoms beyond the second nearest neighbors 
distance and includes interaction of Co with first and second Pt nearest neighbors. 
Herein, we assume that the given supercell is sufficient for Co Ml calculations, which 
is mainly a local quantity. 

We use the LSDA+U method implemented in the full-potential linearized 
augmented plane-wave (FP-LAPW) method including spin-orbit coupling (SOC) 
[9] [16]. When SOC is taken into account, the spin is no longer a good quantum 
number, and the LSDA+U total-energy functional contains additional spin-off- 
diagonal elements of the on-site occupation matrix n miaijm2(T! [6J. The LSDA 
contributions to the effective potential (and corresponding terms in the total energy) 
are corrected to exclude the non-spherical interaction. It helps to avoid the ci-states 
non-spherical Coulomb and exchange energy "double counting" of d-states in LSDA 
and "+U" parts of the effective potential and also corrects the non-spherical self- 
interaction of the d-states. 

In the self-consistent calculations we used 48 special fc-points in combination 
with a Gaussian smearing for the fc-point weighting. A quasi-2D Brillioun zone 
(BZ) with k z = was adopted in order to simulate the 2D-character of the 
problem, notwithstanding that the supercell calculations themselves are inherently 
three dimensional. The "muffin-tin" radii used are Rmt =2.2 a.u. for Co and 2.5 
a.u. for Pt and R'%° T X K max —7.7, with K max the cut-off for the LAPW basis. 
The Coulomb-[/=2 eV and exchange- J=0. 9 eV were chosen which are in the range of 
commonly accepted values for 3d-metals. In principle, U can be calculated by linear- 
response LSDA procedure [11] or from GW [7], both yielding the values ~ 2 eV. As for 
exchange- J, it is not affected by solid-state screening and equal to LSDA calculated 
Stoner exchange parameter. The spin quantization axis is fixed along the out-of-plane 
z-direction. 

First let us make a comparison with the results of KKR for unrelaxed geometry 
(see Table I.). The LSDA calculated values of Ms and Ml agree quite well. However, 
the Ml per d-hole is somewhat bigger in KKR than in FP-LAPW calculations since 
d-shell occupation is bigger in KKR (7.6) than in FP-LAPW (7.2). It is probably 
due to the difference in the radius of ASA-spheres used in KKR and MT-radius in 
FP-LAPW. Since in KKR calculations [T] they use the same radius for ASA-sphere for 
the "big" Pt-atom and the "smaller" Co-atom, this difference can become significant 
and affect both the charge and spin density distributions. 

The XMCD experiments measure not the Ml itself but the Ml per d-hole ratio 
Ml/u^. The Ml/u^ calculated in KKR and FP-LAPW for unrelaxed geometry and 
making use of LSDA is a factor of two smaller than one measured by XMCD. It was 
already mentioned above that Ref. [1] proposed the use of the Brooks OP with reduced 
Br in order to improve the agreement with experimental data. Here we show that the 
Br reduction alone does not solve the problem. Rather it attempts to compensate for 
limitations of the calculations without taking account of structural relaxation in open 
systems including an important class of 3d-adatoms and clusters. 

Next, we turn to the salient aspect of our investigation, the LSDA+U calculations. 
When the Co-Pt inter-atomic distance as for pure Pt is used, the calculated Ml and 
Ml/^ are quite big (see Table I.). By varying dco-Pt and minimizing the total 
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energy we find the equilibrium dc ~pt ~ 3.48 a.u., i.e. reduced by almost 20 %. Since 
we do not perform a full relaxation and the amount of Pt in our supercell is relatively 
small, we can not claim that this will be the correctly optimized dc -Pt- Nevertheless 
it is reasonable to assume that calculated dco-Pt distance is approximately correct 

The change in dco-Pt has strong effect on Ml and Mhjnh (shown in Table I.). 
The Ml / nu becomes fairly close to the experimental value and the agreement for Ml 
is also substantially improved. We should keep in mind that the "experimental" value 
is given as a product of measured M^/n^-ratio times the KKR-calculated of 2.4. 
Making use of the LSDA+U calculated 71^=2.9, we obtain the "experimental" Ml of 
1.34 ± 0.12 /ib, which is in good agreement with the LSDA+U calculated value. 

To understand how the enlargement of the Co Ml in the LSDA+C/ approach 
comes about we consider the spin and orbitally resolved 3d densities of states (tiDOS), 
which are shown in Fig. 2. The spin-resolved eflDOS (see Fig. 2(a)) reveals a substantial 
narrowing of the band width from ^6 eV for hep Co to ~4 eV for the Co ad-atom 
as well as a moderate increase in the spin-splitting, as is expected for the reduced 
Co coordination. The spin-down DOS is split at the vicinity of Ep. When dDOS is 
resolved in terms of cubic harmonics (see Fig. 2(b)), it becomes clear that the spin- 
down peak below Ep posses e g : 3z 2 — r 2 -orbital character while the spin minority 
d-holes are of e g : x 2 — y 2 and t 2£ rOrbital character. 

Since the spin-up Co <i-band is fully occupied, only changes of the spin-down 
band are essential for the Ml enhancement. The m;-resolvcd Co-cQDOS is shown in 
Fig. 2(c) for the spin-up and spin-down channels. The major contribution to the 
increase of Ml originates from \m s — — i;m; = +2) orbital. The Ml enhancement is 
brought about by in-plane spin-down x 2 — y 2 and xy orbitals and much less affected 
by out-of-plane xz,yz orbitals. The spin-down 3z 2 — r 2 (~ |to/ = 0)) orbital does not 
contribute to Ml- This out-of-plane 3z 2 — r 2 orbital is the most localized due to the 
smallest overlap between Co-3<i and Pt-5d electrons. 

It is necessary to mention that our analysis can not be regarded as truly ab-initio 
due to the use of external Coulomb-?/. Herein we make use of a "commonly used" 
value of U=2 eV while it can be at least in principle obtained from constrained LSDA 
calculations [11]. With increase of U the Ml value will increase, and with decrease of 
U it will decrease. Nevertheless our results show quantitatively the role of Coulomb-?/ 
in the Ml formation. 

Also, we did not consider here the magneto-crystalline anisotropy (MAE) induced 
by Co-adatom. In contrast to Ml which is mostly a local property of the Co atom, 
the MAE will consist of contributions from the Co atom as well as the Pt neighbors 
due to strong Pt atom SOC [8]. Most probably, the quantitative studies of the MAE 
will require a bigger supercell; this is the subject of further work. 

Still we can make a rough estimate for the contribution of the Co adatom into 
the MAE. When spin is rotated from the z-axis (out-of-plane) to the x-axis (in-plane) 
direction there is only a little change in the value of the Co atom Ms, from 2.14 /j,b 
(z-axis) to 2.16 /is (x-axis). The change in Ml is substantially greater, from 1.58 hb 
(z-axis) to 1.42 /j,b (x-axis). Indeed, this strong anisotropy in Ml paves the way for the 
strong MAE. Qualitatively, the Co atom contribution to the MAE can be estimated 
making use of Bruno's relation [13] MAE[= (E x - E z )] « -£/4(Mj? - M|), where £ is 
the SOC constant (76 meV for Co-adatom). In the LSDA+U calculations, we obtain 
the MAE of ~ 3.2 meV/Co which is somewhat smaller than the experimental value of 
9.3 ±1.6 meV. A similar estimate for the LSDA calculations yields the MAE of 2.0 
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meV/Co. While our estimate gives the MAE which is smaller than the experimental 
data, it is exceptionally large compared with other Co-based materials: a few tenth 
of meV for Co/Pt and Co/Au multilayers, and 2.0 meV for Co monatomic wire [13]. 

We note that the orbital moment enhancement has been recently investigated 
in Ref. [14] for the case of the Co monatomic wire on the Pt(lll) surface step edge. 
Making use of XMCD the experimental value of the Co atom Ml = 0.68 ±0.05 /is was 
found , which is somewhat smaller than for the Co adatom case. Also, it was shown 
that LSDA yields the Co monatomic wire Ml ~ 0.16 which is substantially smaller 
than the XMCD experimental value (see e.g. Ref. [H]). An account of Coulomb-C/ is 
increasing the Ml value to 0.45 [1b |16j . improving substantially the agreement with 
the experimental data. . 

In conclusion, employing correlated band theory LSDA+U calculations we have 
provided a microscopic picture of the anomalous enhancement of the Co-adatom 
orbital moment. It is found that two major effects need to be included in order to 
essentially improve the Co orbital moment: (i) a correct LSDA+U orbital polarization 
due to the Coulomb-?/ and (ii) structural relaxation of the Co-Pt interatomic distance. 
The calculated value of Ml is found in fairly good agreement with experimental 
XMCD data [1] when those effects are taken into account. 
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Table 1. Spin (M s ), Orbital (Mi) magnetic moments (in /ig), and Orbital 
moment per d-band hole (n^) for a Co atom on Pt(lll) resulting from the LSDA 
and LSDA+U calculations, 



Co ad-atom/Pt(lll) KKR [1 


M s 


Mi 


Mi/n h 


LSDA 




2.14 


0.60 


0.25 


LSDA+OP/2 




2.14 


1.50 


0.63 


CoPtia 


FP-LAPW 


M s 


Mi 


Mt/nh 


LSDA 




2.18 


0.57 


0.20 


LSDA+U 


Uco = 2.0 eV, J Co =0.9 eV 






unrclaxed 


dco-pt=^- 27 a.u. 


2.23 


2.07 


0.70 


relaxed 


dco-pt=3.48 a.u. 


2.14 


1.58 


0.54 


Experiment 


XMCD [1J n h 


= 2.4 


1.1±0.1 


0.46±0.04 


Experiment 


XMCD (n h = 


2.92) 


1.34±0.12 






Figure 1. Schematic crystal structure of a model super-cell: (left) general view, 
(right) top-view with the Co-adatom in the fee position. 




Figure 2. DOS for a Co adatom on Pt(lll): (a, top) Spin-resolved Co-atom 
c(DOS; (b,left) Co-atom cJDOS resolved in cubic harmonics; (c, right) Co-atom 
ciDOS resolved in complex harmonics. 



